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Abstract 

We consider branching of baby universes off parent one in (l+l)-dimensional 
dilaton gravity with 24 types of conformal matter fields. This theory is equiv- 
alent to string theory in a certain background in D = 26-dimensional target 
space, so this process may be also viewed as the emission of a light string state 
by a heavy string. We find that bare energy is not conserved in (1+1) dimen- 
sions due to the emission of baby universes, and that the probability of this 
process is finite even for local distribution of matter in the parent universe. 
We present a scenario suggesting that the non-conservation of bare energy may 
be consistent with the locality of the baby universe emission process in (1+1) 
dimensions and the existence of the long ranged dilaton field whose source is 
bare energy. This scenario involves the generation of longitudinal gravitational 
waves in the parent universe. 
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1 Introduction 



Generally covariant (l+l)-dimensional theories provide convenient framework for con- 
sidering various suspected properties of quantum gravity in (3+1) dimensions (for 
reviews see, e.g., refs.jl], 0]). In particular, the long standing issue of the possible role 
of topology changing transitions and baby universes 0, £|, || ^ [7|, |8j may be naturally 
discussed in this framework. A special feature of (1+1) dimensions which is useful for 
the study of baby universes/wormholes is that some models admit their interpretation 
as string theories in higher dimensional target space. 

The simplest model of this sort is literally the theory of closed strings in the 
Minkowski target space of critical dimensions 0, [II], |TTJ, jl^fl . Indeed, macroscopic 
and microscopic string states may be interpreted as (l+l)-dimensional parent and 
baby universes, respectively. One feature inherent in that model is that the emission 
of a baby universe always requires non-conservation of bare energy in the parent 
universeQIlT], |16j . This non-conservation does not lead to any drastic consequences in 
the simplest stringy model; in particular, the rate of emission of baby strings is finite 
in the limit of infinite size of the parent string |T2|. 

However, one important feature present in (3+l)-dimensional gravity is missing 
in the simplest stringy model. Namely, in (3+1) dimensions there exists long ranged 
gravitational field whose source in energy-momentum (Newton's gravity law), while 
there is no such field in that stringy model. Intuitively, one may suspect that the 
existence of the long ranged field associated with energy and momentum may be 
an obstacle to energy non-conservation in local processes like the emission of baby 
universes. To address this issue, more refined model than that of closed strings in 
critical dimensions, is needed. 

A particularly simple (l+l)-dimensional model where the mass (energy-momentum 
of matter fields) produces long range effects, is the dilaton gravity with matter that 
has been widely discussed from the point of view of black hole physics 0] (for care- 
ful analysis of the notion of ADM mass in that model see refs. |18|, [19||). Here we 
take a different attitude and consider the emission of baby universes, so we simplify 
the model as much as possible. In particular, we set the number of matter fields 

1 It has been argued that the emission of baby strings should lead to the loss of quantum 
coherence for one-dimensional observer at the parent string. Independently, it has been argued 
on general grounds jl3) that the energy non-conservation is inevitable in modifications of quantum 
mechanics allowing for the loss of quantum coherence (see, however, refs. (l4|, ^5)). So, energy 
non-conservation in the stringy model of ref. may not be too surprising. 
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equal to 24 and the (l+l)-dimensional cosmological constant to zero. As stressed 
in refs. p0|, 21] , this model is equivalent to bosonic string theory (in critical dimen- 
sion D = 26 of target space) in linear Dilaton background (to distinguish between 
dilaton fields in (l+l)-dimensional world and in D-dimensional target space, we call 
the former "dilaton" and the latter "Dilaton", respectively). Hence, the emission of 
a (l+l)-dimensional baby universe by a parent universe has an interpretation from 
the D- dimensional point of view as the emission of a light string state by a highly 
excited string state, in complete analogy to ref. [[□]]. This process, in the leading or- 
der of string perturbation theory, is tractable both qualitatively and quantitatively; 
in particular, one is able to analyze whether it is accompanied by non-conservation 
of (bare) energy in (1+1) dimensions and whether its rate is large (unsuppressed) 
when baby universes are emitted locally in (l+l)-dimensional parent universe. The 
discussion of these points is the main purpose of this paper. 

The outline of the paper is as follows. In section 2 we describe the model and 
some of its classical solutions in (1+1) dimensions. For technical reasons, the quan- 
tum version is conventiently constructed for the case of closed (l+l)-dimensional 
universe, so we present in section 2 some classical solutions in the closed universe. 
This discussion will be useful to understand that the (bare) energy of localized dis- 
tributions of matter fields still produces long ranged effects, at least in some gauges, 
even though the total energy of the closed universe is always zero. In section 3 we 
outline the quantum version of this model, which is known for some time (see, e.g., 



refs. [p^, [23"| , |24] , pqj), construct the states of parent universes (DDF states) and ver- 
tex operators corresponding to the emission of baby universes. Section 4 contains 
the main results of this paper. We consider the simplest DDF state of the parent 
universe, which can be interpreted as containing just two dressed matter "particles" , 
and analyze the emission of baby universes by this state in the lowest order of string 
perturbation theory. We find that this emission always occurs with non-conservation 
of energy of matter in (l+l)-dimensional parent universe. We then proceed to the 
explicit calculation of the emission rate. Surprisingly, we find that the rate is unsup- 
pressed even for localized distribution of matter in the parent universe, in spite of 
the long range field this matter produces. In section 5 we conclude by presenting a 
scenario showing that the non- conservation of bare energy of matter may be consis- 
tent both with locality of the emission process and with the presence of long ranged 
field; in this scenario, the ADM mass is conserved at the expense of the generation 
of "longitudinal" gravitational waves due to the emission of a baby universe. 
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2 Model and classical solutions 



2.1 The model 

The action for the simplest version of (l+l)-dimensional dilaton gravity with confor- 



mal matter can be written in the form similar to ref . P6 



S = ~J d 2 a^ (~<f>R + 9 aP d a f%f^j (1) 

where <fi is the dilaton field, f l are matter fields, i = 1, . . . , 24, and 7 is a positive 
coupling constant analogous to the Planck mass of (3+l)-dimensional gravity. The 
coupling constant 7 may be absorbed into the dilaton field, but we will not do this for 
book-keeping purposes. Both in infinite space and in the closed (1+1) -dimensional 
universe the field equations are simplified in the conformal gauge 

g a /3 = e 2p Vap 

where r\ is the Minkowskian metrics in (1+1) dimensions. In this gauge, the fields p, 
cj) and /* obey massless free field equations. There are also constraints 

- ~ 7 2 (d±<j>d ±P - \dl^ + ~ (d ± f) 2 = (2) 
ensuring that the total energy-momentum tensor vanishes. 



2.2 Solutions in infinite space 

Let us outline some classical solutions in this model. We begin with the case of 
infinite one-dimensional space, a 1 6 (—00, +00), and consider localized distributions 
of matter. In this case one can further specify the gauge and choose 

P = (3) 

so that the space-time is flat. Equation (Q) then determines the dilaton field for a 
given matter distribution. Indeed, the solution to eq.(0) is, up to an arbitrary linear 
function of coordinates, 

0(a) = 0+( ( x + )+0_(a_) 

with 

*± = ~f da' ± \a ± - a' ± \ (d ± f) 2 (a' ± ) (4) 
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Hence, the energy-momentum of matter fields produces long ranged dilaton field 
which has linear behavior at large \a |. In particular, the ADM mass can be defined 
as follows, 



Padm — — — 
In 



da 1 



cr 



dcj) 



rV 



-oo 



(5) 



In virtue of eq.(P|) it is equal to 



Padm = da e M (cr) 

J — oo 



where 



1 

2^ 



(<9 f) 2 + (dj) 2 

is the energy density of matter which we will often call bare energy density 0. 

The dilaton field produced by two lumps of matter of equal energy and opposite 
momenta, moving towards each other (or from each other) with the speed of light, 
is shown in fig. 1. Needless to say, the linear dependence of <p on a 1 at large | cr 1 1 is 
nothing but the Coulomb behavior of long ranged field in one-dimensional space. In 
this respect the dilaton field in (1+1) dimensions is analogous to gravitational field 
of Newton's law in (3+1) dimensions. 



2.3 Solutions in closed space 

For technical reasons, the quantum version of this model is conveniently formulated 
in closed one-dimensional space. So, it is instructive to consider classical solutions in 
the closed universe. Let us study the classical theory (|TJ) on a circle 

The absolute length of the circle is irrelevant as we are dealing with scale-invariant 
action; what will matter is the relative size of matter distribution to the length of 
the universe. The gauge (0) is no longer possible in the closed universe; the closest 
analog is the gauge in which the universe contracts (or expands) homogeneously, 

p = --P(-V° (6) 

7 

where P+) is some constant; our choice of normalization and notation will become 
clear later. In this gauge, the constraints (|2|) may again be used to determine the 

2 Note somewhat unconventional factor 1/tt in the matter action in eq.(0). 
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dilaton field for a given distribution of matter, provided the total spatial momentum 
of matter vanishes, 



/■+7r/ 2 /-+7r/ 2 

/ da 1 (<9 + f) 2 = / da 1 (<9_f) 

J-n/2 J-n/2 



tt/2 
/2 



(7) 



The constraints (0) in the gauge (^) read 



-i 7 p(-)5 ± 0-i 7 2 40 = i(9 ± f) 2 (8) 

Equation (|7|) is an immediate consequence of these constraints and the periodicity of 
the dilaton field in a 1 . 

For a given matter distribution, the solution to eq.(^) which is periodic in a 1 with 
period n is 

(9) 
(10) 



(o°,(7 1 )=6 + (a + ) + 6Ja.] 



where 



£t(<r±) = / da ± G(a ± ,a' ± ) l -{d ± f) 2 
and the Green function of eq.(Q) obeys the periodicity condition 

d u G{u + it, u) = d u G(u, u) 

At u G (— 7r/2, 7r/2), «' G (— 7r/2,7r/2) we have explicitly 



«') 



7PH sinh(7rP(-)/27) 



exp 



7T 



p(-) 



20 



-el tt — u 



exp 



pH 
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[u — u 



- 1 



where e(u — u') is the usual step function. 

It is instructive to consider the case 7 > 1 (large "Planck mass") and study two 
narrow pulses of matter moving left and right and colliding at cr 1 = 0. These pulses 
may be approximated by the delta-function distribution, 

7T 



(d ± f) 2 = -»5(a ± ) 



(12) 



where the normalization is such that the constant fi would coincide with the ADM 
mass had the universe infinite size, 



da 1 — 
2tt 



(<9 f) 2 + (d 1 ff 



In this case eq.(|10|) has particularly simple form (again at |cx±| < 7r/2) 



p(+) 
'~27~ 



7T 



2 7 



2 7 2 



(13) 



(14) 



where 

P' +> = ^ (15) 
When the pulses are close to each other, one has for the total dilaton field at |cr±| <C 1 

P(+) 7T 

-o°- 7 ^ti\°+\ + W-\) + - (16) 
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The first term in eq.(|T6D describes spatially homogeneous component of the dilaton 
field, while the second term shows precisely the same Coulomb behavior as in the case 
of infinite space. In fact, the latter term coincides with the expression (Hj) (for narrow 
matter pulses as defined in eq.(|T2"D). The terms omitted in eq.(|T6|) are of order of 
7~ 2 and they become important at |er±| ~ 1; in particular, they ensure that 4>(cr°, cr 1 ) 
flattens out and has vanishing spatial derivatives at a 1 = ±7r/2. The behavior of the 
dilaton field generated by narrow pulses of matter, which are close to each other, is 
schematically shown in fig. 2. 

Even though the ADM mass is, strictly speaking, zero in the closed universe, 
narrow pulses of matter produce the dilaton field that shows Coulomb behavior not 
far away from the pulses. In this sense one can still use the notion of bare energy 
in the gauge (](]). This bare energy (the energy of matter) can be observed by a 
one-dimensional observer by measuring the dilaton field outside the pulses but at 
distances small compared to the size of the universe. With these reservations, the 
formula for the ADM mass, eq.(|||), still makes sense (when all events occur and fields 
are studied in a small part of the universe). Clearly, all these observations apply only 
to those matter distributions whose conformal size, a 1 ~ r puise^ j g sma }} compared 
to 7r, the conformal size of the universe. In other words, we will be interested in 
considering large wave numbers, 

n -j- > 1 (17) 

repulse v ' 

The gauge © is not the only useful one in the closed universe. At large 7, one 
can choose the gauge 

p(+) 

V° (18) 



7 

where prime is used to denote the quantities in this gauge. In general, the coordinates 
in the two gauges are not too different, 

4 = ^ + 0(7^) (19) 
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The coordinate transformation has a particularly simple form in the case of two 
narrow pulses; it follows from eqs. flllD and ( |l5f) that in that case 

<4 = *± + ^) k±l " ^T)^ + °(7~ 2 ) (20) 
Furthermore, for two narrow pulses one obtains (at \a±\ < 7r/2) 

This expression, of course, solves the constraint (§) for matter distribution 

(* ± f)V ± ) = |/i*(<4) (22) 

The fact that the matter distributions (O) and ( P2"[ ) are essentially the same in the 
two gauges, is an immediate consequence of eq. (|19D . 

If the pulses of matter are not infinitely narrow, the above features remain valid 
qualitatively. The matter distributions in the two gauges (^) and flIB|) are the same 
up to corrections of order 7 _1 ; in particular, the spatial sizes of the lumps differ only 
by a factor (1 + 0(7 -1 )). The scale factor p'(a') changes rapidly in the regions of 
non- vanishing energy-momentum density of matter, and has the form shown in fig. 3, 
the depth of the well being proportional to the total energy of mater. We will have 
to say more about this gauge in section 5. 



3 Quantum states and vertex operators 

The model is quantized exactly in the same way as bosonic string theory in Minkowskian 
target space of D = 26 dimensions. One introduces the notation 

<f>=--(X° + X 1 ) 
7 

p = -i(X°-X 1 ) 

7 

/\...,/ 24 = X 2 ,...,X 25 (23) 

Then X m (<t), li = 0, . . . , 25 are canonically normalized free two-dimensional fields and 
the classical constraints @ become 

\d ± X,d ± X» - ^ 7 <9|X« = (24) 
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where the summation is performed with Minkowskian D-dimensional metrics, r]^ u = 
(— 1, +1, . - - , +1) and 

X& = X°±X 1 

Upon quantization, the left- and right-moving components of these fields are decom- 
posed in the usual way, 

= \^ + \P>°- + \ E 

where and a£ are the standard oscillator operators with string normalization. 
According to eq. (|24]) , the Virasoro operators are 

1 



;^ 2 + E 

fc>0 



a -k a k 



L m = + \ E : «-k ■ +^7^"L +) (25) 

and similarly for L and L m , where affi = ± a^. Equation (p5|) shows that the 
model is equivalent to the string with background charge [^, ^3], ^4], ^5] or, in other 



words, to the bosonic string in linear Dilaton background [20, 21 



= 7(e^) = 7x (+) (26) 

with e M = (—1, 1, 0, . . . , 0) being a light-like vector in D-dimensional target space. 

It is clear from eq. ( ^5l) that the spectrum of states in the target space is the same 
as the spectrum of the bosonic string in trivial background, 

M 2 = 8(n-1), n = 0,l,... 

The vertex operators are, however, slightly different. For example, the tachyon vertex 
operator of conformal dimension (1,1) is 

V(Q) =: e ^ x ^ x(+) : (27) 

with Q 2 = 8. This modification can be understood as being due to the Dilaton 
background fl26|) . Indeed, the effective action of tachyon field T in flat target space- 
time and in linear Dilaton background $(x) is 

1 

2 



£> „-2$ 



d x e 



; (<9 M T) 2 + 4T 2 + cT 3 + 0(f 4 
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By introducing the field T = e *T one rewrites this action in the following form, 

d D x 



-l(^T) 2 - ^T 2 + ce^T 3 + 0(T 4 



where M| = — 8 + (c^) 2 , i.e., M| = —8 for light-like e M in eq . (|26|) . In this notation 
the kinetic and mass terms are conventional, while the trilinear vertex is propor- 
tional to e*, precisely as required by eq. (|27D . This argument can be generalized to 
interactions of the fields other than the tachyonpO], 



Let us now consider highly excited string states (parent universes). They are 
conveniently constructed by making use of the DDF operator sp8|. In the light-like 
Dilaton background, the simplest choice of the DDF operators is 



-7T/2 



a: 



A', 



-7T/2 
+7T/2 



7T 



7T 



exp 



exp 



4in 



-7T/2 



d + Xl(a, 



(2? 



with the same light-like vector e^. Here i = 2,..., 25. These operators obey the 
usual oscillator commutational relations, and their commutational relations with the 



(29) 



Virasoro operators [|11[] ensure that the state of the form 

A\ 1 ...AX 3 -A^ i ...A^ t \V) 
is a physical state provided that \V) is the physical tachyon state and 

m + . . . + n s = hi + . . . + h t (30) 
The DDF operators (E3) are similar to those introduced in ref.pSfl. Another choice 



of DDF operators in this model has been considered in refs.[23, 25 



To make contact with the classical analysis of section 2, consider suitably modified 
coherent states 

l*P> = P L =L Q ( E \f l nX-n + E XP) (31) 

\n>0 ' l n>0 " / 

Here /* and /| are c-number amplitudes and P Lq= i is a projector onto the subspace 
of vectors obeying 

(£o-L )|*) = 
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This projector is needed to ensure the validity of eq.(pO[) term by term in the expansion 
of \^ P ). 

Let us impose the condition that 

E^ = EK ( 32 ) 

n h 

and take the amplitudes f n and to be large. Consider now matrix elements of the 
form 

(y pl \0 M ({a + };{cT-})\*p) (33) 

and also 

(* P ,\d-Xg\<j'_)0 M ({a+y, W-})\y P ) (34) 

(* P ,\d + X^(a > + )0 M ({a + } ] {<r_})|¥p> (35) 

where the operators Om are products of matter fields, in general at different points 
in (l+l)-dimensional space-time, 

M = d+Xiiai) . . . d + X^{a\) ■ d-X*(oi) . . . d_X%(al) (36) 

and, as before, we take ij, . . . ,j g — 2, . . . , 25, so that the operators Om indeed contain 
matter fields only. The matrix elements (|33|) are then the correlators of matter fields 
in the coherent state (|3l[) , while the matrix elements (34) and (|35|) are the correlators 
of dilaton and metric fields with matter. 

These matrix elements are calculated in Appendix with the following result. Up 
to small corrections and trivial normalization factor, they coincide with classical cor- 
relators 

f —o M ({a ++ ey,w--e})^ / ^^(Km^+o) (37) 

J-tt/2 7T J-tt/2 7T 



and 

— d + Xt {± \a + + e)0 M ({a + + W- ~ ?}) (38) 

J-tt/2 IX 

(and similarly for eq.(34)) respectively, where M is given by eq . (|36D with classical 
matter fields 

d + X d >\o + ) = \P 1 + £ (f>' 2ina+ + fn* 2 ™" 

1 n>0 

The classical field X cl ^ is defined by 

Q X cZ,(+) = pW 
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while the field d±X cl,<y ~' > is to be found from the classical constraints (|24]). In short, 
the matrix elements like eq.(^) and eg. (p5|) are equal to the corresponding classical 
expressions in the gauge (IS), integrated over translations in one- dimensional space. 



In this way the classical picture is restored; the particular choice of the DDF operators, 
eq. (|28|) , corresponds to the gauge choice (|I~8"D. Note that eq.(p2[) is precisely the 
classical constraint Lq = Lq written in this gauge. 

Hence, the DDF operators (|28| ) correspond to creation and annihilation of dressed 
matter excitations in (l+l)-dimensional space-time in the gauge <fi = const ■ a . The 
corresponding wave numbers are n and n, respectively. In what follows we will call 
(somewhat loosely) these excitations as "dressed particles" in the (l+l)-dimensional 
universe. 



4 Non-conservation of bare energy in (1 + 1) di- 
mensions and emission probability 

4.1 State of the parent universe 

In this section we consider the simplest DDF state 

\%n,i,j) = J ^{P)\P,n,i,j) (39) 

where 

\P,n,i,j) = -A i _ n AL n \ , P) (40) 
n 

and ^{P) is the wave function of the center-of-mass motion in target space, in mo- 
mentum representation. The normalization factor 1/n in eq. (f4"0|) is chosen in such a 
way that the state \P,n,i,j) has the usual D- dimensional normalization (recall the 
string normalization of the oscillator operators, [A^, = nS^). According to the 
discussion in section 3, we interprete this state as the state of a parent universe with 
two dressed matter particles (one left-moving and one right-moving) of equal wave 
numbers n, in the gauge ([18]). The normalization convention in eq.(|4T)|) corresponds to 
"two particles in entire one-dimensional space" normalization of (l+l)-dimensional 
quantum field theory. These particles have equal bare energies and opposite bare 
momenta in (1 + 1) dimensions, 

Cleft = bright = ^ (41) 
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Pright = -Pleft = 

We are interested in the limit (see eg. (|17[)) 



n — ► oo 

By making superpositions of the states fl39|) with different wave numbers n one can 
construct states with localized distributions of matter in one-dimensional universe. 
This generalization is straightforward, so we stick to the state with fixed n. The 
necessity to consider wave packets (^) in target space and not just plane waves ( f4"0| ) 
is due to the Dilaton background that increases indefinitely as — > oo: amplitudes 
of processes involving plane waves would be divergent in this background. 

It is a matter of simple algebra to see that the total D- dimensional momentum of 
the state (gDJ) is 

P k = v k^ = 2, . . . , 25 
p(+) = pO + pi _ 

P(-) = po _ pi = p(-) _ (42) 



Hence, from Z)-dimensional point of view, the state (HO) is interpreted as an excited 
string state at the n-th level with the mass 

Ml = 8n - 8 (43) 

(recall that \V) is the tachyon state). We consider for definiteness this state in the 
center-of-mass frame, 

p fc = 0, A; = 2, ...,25 
p(+) = p(+) = Mn 

P { - ] = M n (44) 

although our discussion can be straightforwardly generalized to other frames. In terms 
of the wave packets (p9|), eq.(|j) nieans that the wave function of the center-of-mass 



motion, \ty(P)) is peaked near the values determined by eq. flSil) . 



4.2 Emission of baby universe: energy non-conservation in 
(1+1) dimensions 

The parent universe in the state (P9|) can emit a baby universe — a universe with no or 
small energy of matter particles. In D-dimensional language this process corresponds 
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to the emission of a low lying string state (tachyon, Dilaton, etc.) into D-dimensional 
target space. 

Let us study whether this process always occurs with non-conservation of bare 
energy in (1+1) dimensions. Since bare energy coincides with the level of the string, 
eq. fl41|) , we are interested in the change of the D-dimensional mass of the highly 
excited string due to the emission of the low lying string state. At first sight, the 
presence of the time-dependent Dilaton background in target space might give rise to 
non-conservation of .D- dimensional energy and momentum in the emission process. 
In that case the emission of a low lying state would not necessarily require the change 
of the level of the heavy string; in other words, the emission of a baby universe 
would not necessarily require non-conservation of bare energy in (1 + 1) dimensions. 
Surprisingly, we will see in a moment that this is not the case: D-dimensional energy 
and momentum are in fact conserved in the presence of linear Dilaton background 
exactly as they do in flat target space with no background. Hence, the emission of a 
baby universe always occurs with non-conservation of (l+l)-dimensional bare energy 
of matter in the parent universe. 

The argument presented below is fairly general; the particular form of the ini- 
tial state, eq. (|39|) , is unimportant. Non-conservation of energy in (l+l)-dimensional 
parent universe due to the emission of baby universes is a generic property of our 
model. 

The argument for the conservation of D-dimensional energy-momentum in the 
presence of the linear Dilaton background is conveniently presented by considering 
a theory of three scalar fields $1, $2; $3 with cubic interaction and exponentially 
changing coupling. Let the cubic coupling be 

g(x) = exp(I>") (45) 

with real constant T^. Let us take $i to be the heaviest (mass Mi), and $2 to have 
the mass Mf which for simplicity is close to Mj (but smaller than Mj). Let the field 
$3 be light. We are interested in the process when a particle $1 emits a particle $3 
and becomes a particle $2- 

Because of the unboundedness of the coupling ((!§), it does not make sense to 
consider plane waves of particles $1 and $2- Instead, one has to use wave packets. 
Namely, consider the amplitude of the decay of a wave packet ^/j(x, t) describing the 
state of the particle $1 into a wave packet f of the particle $2 plus a particle $3 
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which has fixed energy-momentum Q 



A = J dxdt g(-K, t)*i(x, t)^}(x, t) exp(-iQ ^x**) 



(46) 



We will see in the next subsection that the amplitudes of string decays have similar 
form. 

Let us specify the form of the wave packets and ^f. Let us consider for 
definiteness the non-relativistic regime and take wave packets narrow in momentum 
representation. Let us furthermore neglect the dispersion of the wave packets with 
timef] (this can be achieved by confining the particles in moving potential wells). 
Thus, we take 



exp — iMit — 



— P? + iPiX - — (x - Vit) 2 

2M, ; 1 2 y 1 



(47) 



where a, the width of the packet in momentum space, is small, and Vj = Pj/M; is the 
velocity of the particle, which is also assumed to be small for simplicity. Similarly, 
the wave packet has the form ([4"7|) with Mj, P«, v, substituted by Mf, P/, v/. 

The integral ( [46]) is then Gaussian, and is straightforward to evaluate. The result, 
up to a pre-exponential factor, is 



.4 



exp 



2a 2 



-*i(AP, AE) + Fi(T, T ) + iF 2 (AP, T, AE, T )) 



(4J 



where 



AE 




are amounts of non-conservation of momentum and energy and 



F!(AP,A£) 



F x (r,r 



0^ 



(AP) : 

4 

(r) 2 



+ 



i 



i 2 



+ 



(v 4 - V 



A£?--(v 4 



V/ )AP 



f) 



3 The dispersion of wave packets with time would complicate the analysis considerably. The 
coupling changes so rapidly that in the case of spreading wave packets the interaction often occurs 
in space-time very far away from the centers of the wave packets. In that region there are mostly 
modes with momenta quite different from the central values Pi or P/. Hence it is difficult to 
separate the effects of finite widths of the wave packets in momentum space from possible effects of 
momentum non-conservation. 
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The explicit form of F2 is not important; it is sufficient to note that both F% and F2 
are real. The only important property of F\ is that it is positive-definite and vanishes 
iff AP = AE = 0. 

The imaginary part of the exponent, iF 2 , in eq.(^Sj) is unimportant and cancels out 
in the probability. Then the probability factorizes into a term depending on AP and 
AE and a term containing T M . The latter term is nothing but the overlap of the initial 
and final wave packets with the weight g(x). More importantly, since Fx(AP, AE) 
is multiplied by the large factor 1/a 2 in the exponent, and because Fi(AP,AE) is 
non-negative and vanishes only at AP = AE = 0, the exponential factor containing 
Fi(AP, AE) ensures conservation of energy and momentum in the limit of small a 
exactly in the same manner as it does in the case of space-time independent coupling. 

The restriction to the non-relativistic case and Gaussian wave packets is, in fact, 
not essential: the same argument goes through for relativistic and non-Gaussian wave 
packets (provided that they do not disperse with time). 

The result that energy and momentum are conserved in spite of space-time depen- 
dence of the coupling is peculiar to the exponential coupling whose exponent linearly 
depends on x M : only in this case the dependences on AP M and factorize. For in- 
stance, if the coupling switches off at infinity in space-time (say, has finite support), 
then the same calculation leads to the usual result that energy and momentum are 
not conserved, and the amplitude is proportional to the Fourier component of the 
coupling, g(AP, AE). 

4.3 Emission rate 

Let us now turn to the actual calculation of the rate of the emission of a baby universe 
(low- lying string state) by a parent universe (excited string) in the state fl39l) . We 
explicitly consider the emission of a tachyon, although the analysis — and results - 
are the same for the emission of a Dilaton or graviton as the baby universe. Let the 
.D-dimensional momentum of the outgoing tachyon be Q^. We first have to specify 
the range of which is of interest for our purposes. 

The emission of a baby universe due to the collision of (l+l)-dimensional particles 
with wave numbers equal to n has a chance to be local if the characteristic conformal 
time of the process of emission, Act , is of order 

A<x°~- 

n 
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This conformal time is related to the /^-dimensional time via 



Ax° = P°Aa° (49) 

The D-dimensional time characteristic to the tachyon emission can be estimated in 
the center-of-mass frame of the decaying string as 

This is the time after which the tachyon is formed and splits off the initial string. 
Hence we are interested in the tachyon energies of order 

go _____ ^ (50) 



where we made use of eqs.(f43l) and fl44]) . Since Q° is large, we neglect the tachyon 



mass where appropriate in what follows. 

We are interested in the process in which the initial string at level n decays, by 
emitting a tachyon, into level n' . According to eqs.(|43D and (|50|), and because of 
energy conservation in D dimensions, we have 

m = n — n'r^n (51) 

We intend to sum up over all final states at given level n' . Usually, this summation 
is conveniently performed by evaluating the forward amplitude shown in fig.4 where 
curved lines denote the tachyon (cf. ref. [0). Because of the space-time dependent 
coupling, the procedure is somewhat tricky in our case. It will be convenient to 
consider first the amplitude with different momenta, as shown in fig. 5, and set P = P' 
and K = Q in the end. The amplitude is given by 

A*{P, P'; K, Q) = k 2 {P\ n, i,j\V(K)AV(Q)\P, n, i,j) (52) 

where k is the string coupling constant, V is the vertex operator (pTj) and A is the 
usual string propagator (recall that the Virasoro operators Lq and Lq coincide with 
conventional ones). Here we consider the plane wave (f40"D as the initial state; the fact 
that we actually have to deal with wave packets like (^|) has been already discussed 
in the previous subsectionf_ 



4 The procedure below involves manipulations with formally divergent integrals, etc. This proce- 
dure can be checked by explicit analysis of amplitudes involving low lying string states only. 
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The evaluation of the amplitude ([52]) is straightforward and parallels that of 
ref.fll2|. One finds 

A^(P,P f ;K,Q) = k 2 J d D xe i{P-Q> e -HP'-K)), e 2 lx M 

x / dzdz {zzf^-^^^^^F'Hz.z) (53) 
J n 2 



where we have kept the integration over the center-of-mass coordinate x and denoted 



with 

B ij {z) 



F i] (z,z) = \i- z \h( Kfl -i-y^)(-Q^'y^)B ij (z)B ij (z) 



(54) 



du du' 1 



2tt 2tt u n+1 (u') n+1 



r(+) 



X (1 -u) 1 - 



nQ( + ) 

U\ P (+) 



f(+) 



n Q (+) 



(l-Ji')7!+T(l-ztf) ^' (+) 



4 V z — u 



1 — it 



^ 1 



+ 5 lj - 



uu 



(55) 



(1 -uu'f 

The integrations here run around small circles in complex u- and w'-planes surround- 
ing the origin. In eqs.(^) and ([)4]) we used the light-like vector e M = (— 1, 1, 0, . . . , 0) 
defined in section 3. 

To extract the decay rate into final states at given level n', we expand F(z, z) in 
a formal series in z, z (we omit superscripts i, j temporarily), 



+oo 

F{z,z)= £ F mm ,z' m z~ m ' 

m,m'=—oo 



(56) 



At m = m' the corresponding integrals in eq.(|54]) have poles, 



dzdz 



irF n 



_|(PM - i 7e M)(-g M - i 7 e M ) - m - 1 

8tT F'min 

(^-Q,-ne,) 2 + (M2-8m) 



(57) 
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These pole terms lead to contributions to the amplitude ([5^) which can be written in 
the following form, 

AV(P,P>;K,Q) = J2^ 2 ^F mm [ d D xd D yd D P f e« p - Q ^ x{+) 

m n J 

c .—iPf(z—y) 



I 2 , 

L n—mJ 



This expression is recognized as the sum of the amplitudes of processes going through 
states with masses M n _ m in a theory with trilinear coupling which exponentially de- 
pends on x^ + \ Hence, the probability of the decay into level n' = n — m is determined 
by F mm . The total probability involves also the overlap between the initial and fi- 
nal wave functions of the center-of-mass motion. This overlap has been considered 
in previous subsection; it leads to the conservation of energy and momentum in D 
dimensions. Therefore, we can set P' = P and K = Q and obtain 

53 \AJ(n ^ n -m;Q)\ 2 = Stt^^F^P' = P\K = Q) (59) 
y J n z 

where A % j denotes the amplitude of the decay into a final state / at level (n — m) and 
a tachyon with momentum Q; the sum in the left hand side runs over all final states 
at the level (n — m). This expression should be integrated over the phase space of 
the two final string states. 

To estimate the integral over the phase space, let us study first the behavior of 
F m m at large n and m. We recall that we are interested in tachyons with energies 
Q° ~ \fn. Let us first consider the generic case, 

Q(+) ^ q(-) ^ Qk ^ ^ (60) 
At K = Q and P' = P we have 

Fmm' = Ft m R m > (61) 

where the form of R 1 ^ follows from eqs. (|54]) , (|55|) and (0), 

r du du' dz 1 
J 2^ 2^2^ u n+1 {u') n+1 z- m+l 

„q(+) 

"(1- 2] (l-w')l P( 

x (1-zX 



(1 - u)(l - zu') 
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X 



u'(l-z) 



+ 8 ij - 



uu 



(l-uu') 2 _ 

In the regime (|60|), i? m can be written as follows, 

R l ^ l = J du du' dz P lj (u,u', z)e~ 



(62) 



where 



S = log u+log u log z — — 

n M, 



n,S(u,u' ,z) 



U , 



(63) 



- log( i — a ) + log(l - -) + log(l - «') - log(l - zu') 



and P w is a pre-exponential factor that depends on n only weakly. Here we made use 
of the fact that p(+) = M n in the center-of-mass frame we consider. 

In the regime flBOD all terms in S are of order one, so the integral in eq.(|63"D can 
be calculated by saddle point technique. At the saddle point one finds 



S = X 



m 



n 



X 



ni 



n 



2Q<+>\ 



M 



(64) 



where 



X {v) = {l-v) log(l - v ) + (l + u) log(l + v) 

Making use of energy-momentum conservation, one obtains (again in the center-of- 
mass frame of the decaying string) 



m (m 2Q( + ) \ m 
n \n M n J n 



so S is always positive. Therefore, we conclude that the decay probability is expo- 
nentially suppressed in the kinematical region (|60|), 



P(n — > n — m; Q) oc e 



-2nS 



Equation (p4|) implies that the decay probability may be unsuppressed in the 
kinematical region different from eq. (|60|) , namely, at 



Q i+) ~ 



(65) 
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i.e., at 



M, 

In this region 



nQ {+) , x 

- 1 (66) 



\nj M n 

Clearly, the saddle point calculation is not valid in this region, so we proceed in a 
different way. To estimate the integral in eq.(62) in the regime fl66f), we make use of 



the following asymptotic formula|12] 



J dz x ... dz k z^ . . . z^ Xak U(l - z p ) a "U(l - z p z q )^U(l - z p /z q y™ 



l\J2 a p+J2 Ppi+Y, "/pi+ k 

oc 



which is valid as A — > oo with a p , a p , (3 pq , 7 Pg fixed. We obtain 

B%<xn 

at large n. We then recall eqs.(|59|) and ( 6l|) and find that at large n 



\A l j(n — > n — m;Q)\ 2 = independent of n 

f 

in the kinematical region (|65|). 

To estimate the total emission probability, we note that in the region (RH) one has 



Q| ~ 1 

where Qt = (0, 0, Q 2 , . . . , Q 25 ). Hence, the probability to decay into the level (n — m) 
is of order 

where E n _ m ~ M n ~ y/n is the energy of the final excited string, and Q° ~ y/n 
according to eq. (|50|) . The number of final states contributing to the decay, Am, is of 
order n, so we have finally 

Pin — > n — m) ~ — == 



This is the decay probability per unit D- dimensional time x°. To obtain its interpre- 
tation in (l+l)-dimensional terms, we recall eq.fl49|) and find that the emission rate 
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per unit conformal time a of (l+l)-dimensional universe is independent of n at large 
n. 

This is our principal result: the rate of the emission of baby universes is unsup- 
pressed at large n, when the emission process should occur locally. This emission rate 
is proportional to the collision rate of two narrow wave packets in one-dimensional 
universe of conformal size it, the proportionality constant being independent of the 
one-dimensional momenta of the "particles" or width of their wave packets and being 
determined by the string coupling constant only. 

5 Discussion and conclusion 

Let us summarize our results for the simplest version of the dilaton gravity with 
conformal matter in (1+1) dimensions. We considered mostly the case of compact 
one-dimensional universe and studied pulses of matter whose size is small compared 
to the size of the universe (i.e., whose wave numbers n are large). At least at the 
classical level these pulses, in the gauge 

p = const ■ a (67) 

produce long ranged dilaton field which is approximately Coulomb at scales small 
compared to the size of the universe. The magnitude of this long ranged field is 
proportional to the energy of matter, which we called bare energy. The notion of 
ADM mass makes sense at these scales and coincides with the ADM mass defined for 
infinite space. 

To construct quantum states, it was convenient to work in a different gauge, 

<f) = const ■ a (68) 

It is important that at large 7, the sizes of matter pulses in the two gauges are 
similar, again at the classical level. In other words, the pulses that are narrow in 
the gauge (p7|) are also narrow in the gauge (|68"D, so the processes we were interested 
in are local in both gauges. Making use of the gauge (|68|), we considered — at the 
quantum level — the simplest state of the parent universe that contains one left- 
moving and one right-moving dressed matter "particles" with large wave number n. 
The collisions of these particles may eventually induce the emission of baby universes. 
If the relevant quantum numbers of the baby universe (D-dimensional momenta 
of the microscopic string) are large enough, the emission process is local in the parent 
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universe. We have found that the emission always occurs with non-conservation of 
energy of matter, and that the probability of this process is finite at large n. 

At first sight, there appears to be a conflict between the locality of the emission 
of a baby universe, and hence the locality of the non-conservation of matter energy, 
and the existence, in the gauge (0), of the long ranged dilaton field whose strength 
is determined by the matter energy. To see that this conflict is only apparent, let 
us present a scenario consistent with both of the above properties. We stress that 
the following consideration is only a scenario, as its confirmation or rejection would 
require the analysis of the final state of the parent universe, which goes well beyond 
the scope of this paper. Also, the discussion below is essentially classical, while the 
actual analysis should necessarily be at the quantum level. 

Let us again consider the collision of two narrow pulses of matter, and choose 
the gauge (|68|). In this gauge the field p and the matter energy density are those 
shown in fig. 3 and fig.6a. If the matter energy was conserved, the final state would be 
characterized by the configuration shown in fig. 6b by dashed lines: the field p between 
the pulses would change from — p to +po where p is determined by the total matter 
energy in the pulses (see eqs.(^Tj) and (fL?)) , flllf )). If the energy is not conserved at 
the moment of the collision (i.e., if the collision of the pulses induces splitting off of 
the baby universe), the height of the matter pulses and, correspondingly, the height 
of the plateau of p are smaller in the final state; this configuration is shown in fig. 6b 
by solid lines. 

Clearly, the process shown in fig.6 may be perfectlly local in (l+l)-dimensional 
space-time. It shows that non-conservation of matter energy does not require non- 
locality. However, this process cannot be transformed into the gauge (|6"T[), as the 
field p does not obey the field equation d a d a p = everywhere in space-time. To 
see what happens if the gauge (|67|) is chosen for the initial state, let us perform the 
gauge transformation that would transform the "conventional" configuraion (i.e., the 
configuration of the conventional process with energy conservation) into the gauge 
(|67D. In the case of infinitely narrow pulses this gauge transformation is the inverse 
of eq. (pUD . Then the initial state is one shown in fig.2 (with p = const -a everywhere), 
while the final state is that shown in fig. 7 by solid lines (only a small region of the 
universe is presented in fig. 7; the final configuration of the conventional proccess with 
energy conservation is again shown by dashed lines for comparison). The final dilaton 
field in this gauge is the same as that of the conventional process; in particular, its 
long range behavior is not affected by energy non-conservation. On the other hand, 
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the field p in the final state is non-trivial and corresponds to longitudinal gravitational 
waves. It is the presence of these longitudinal waves that ensures the validity of the 
constraints after the collision, even though the energy of matter is not conserved 
and the dilaton field does not change asymptotically. Of course, the longitudinal 
p-wave may be gauged away, but the corresponding gauge transformation would be 
non-trivial far away from the collision region, and would also induce longitudinal 
gravitational wave in the initial state. 

In infinite space, the gauge (|68|) cannot be imposed, so we cannot use the ar- 
guments based on fig. 6. However, the final states like those shown in fig.7 are still 
possible in the gauge where p = initially. For these final states to appear, the 
field equations should be violated only in a small region of space-time (where the 
two pulses collide), and the entire process may occur locally. The ADM mass viewed 
from infinite distance is conserved, but this conservation is due to the appearance of 
the longitudinal gravitational waves that compensate for non-conservation of matter 
energy. 

It remains to be understood what part, if any, of the discussion of this paper 
may be relevant to (3+l)-dimensional theories. There exist semiclassical arguments, 
based on the study of fluctuations about [59] and analytical continuation of|3(J the 
Euclidean wormhole solution of ref.||, favouring the interpretation of the wormhole 
as describing the process in which a baby universe branches off and then "flies away" 
in (mini)superspace. This process may be very similar to the one discussed in this pa- 
per in (l+l)-dimensional context. On the other hand, the possible non-conservation 
of bare energy was not explicit in the semiclassical treatment of (3+l)-dimensional 
Euclidean wormholes. An independent problem which can possibly be treated "phe- 
nomenologically" , as we did in this section, is whether non-conservation of bare energy 
in (3+1) dimensions is consistent with locality, and, in particular, whether locality 
requires the generation of longitudinal gravitational waves. We hope our study of 
(l+l)-dimensional toy model will be helpful to understand these problems. 

The author is indebted to A. A. Tseytlin for very helpful correspondence and to 
P.G. Tinyakov for useful discussions. This work is supported in part by the Russian 
Foundation for Basic Research, project 96-02-17449a, by INTAS grant 93-1630-ext 
and CRDF grant 649. 
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Let us outline the calculation of the matrix elements fl33|) and (|35|) in the leading 
order in f n , f„. First, we make the projection onto the subspace L = L explicit by 
writing the coherent state (|5T|) in the following form, 



I*p> 



i/2 rff / 1 1 „ „ \ 

— exp X -fn(OA-„ + E 4fn(^)A-n |P> 



(69) 



where 

fn(0 = f«e 2m « 

fn(0 = he~ 2 ^ (70) 

The representation ( |69"D coincides with eq. ([Jl]) up to normalization. The fact that the 
state (|6~9|) obeys the constraint (L — Lo)|^p) = follows from the commutational 
relations of the DDF operators with L and L , 

n 



[L ,Al 



-A r 



n 



[L Q ,A* n ] = -A\ 

and similarly for A 3 A . 

Consider now the norm of these states. One has 



7T- 



exp (£fj;;e 2m ^) + £f.f* e - 



At large f n and this is a saddle point integral. Taking into account eq.(|32"D we find 
that the integrand does not depend on (£i + £ 2 ), while the saddle point in (£i — £ 2 ) is 
at 

6-6 = (7i) 

Hence we obtain the usual result 



up to a pre-exponential factor. 

Let us turn to the matrix elements (P3|) involving matter fields only and consider 
explicitly the left-moving sector. The DDF operators can be written as follows, 



At 



V2 da. 



x 



(+) 



X 



tt/2 7T 
1 



exp 2m— —t- + Una 



p(+) 



r ; ■ E" 



t A -2j<jo-+ 



(72) 
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Since P l and pW commute with A % nl one can set 

pi _ pi p(+) _ <p(+) 

in the operator O m for calculating the matrix element (|3~3"D. Furthermore, com- 
mute with Om and with all factors in eq. (|72|) , so one can set them equal to zero and 
write effectively 

We have to calculate the matrix element 

^ pl \a\ i ...a\ a -al\...a^ P ) (73) 



with ri, . . . , r s ,pi, . . . , p t > 0, which is a building block of eq. (|33|) . Note that the 
operator ordering in eq.([73|) is in fact not essential at large f n , as is usual in the 
classical limit. One finds for this matrix element 



/ 



{V | exp 



TV 2 



x (+) 

2i(n + ... + r s -p 1 - Vt)p^ 



(RM)|P) 



x/;(6) • • • £(60 • ■ ■ ■ Uti) ■ (RM) (74) 



where we omitted the superscripts ii, . . . ,jt', (RM) denotes the corresponding factors 
due to right-moving modes. All dependence on (£i + £2) in this integral comes from 
the exponents in /*(^) an d f P (£,i), see eq. Q7DD , and similar exponents for right-moving 
components. This makes the matrix element in eq.([74]) equal to (V'\V). The integral 
over (£1 — £2) is still of saddle point structure with the saddle point (|71~P . Hence the 
expression ( [74] ) simplifies and becomes equal to 

(* P ,\* P ). J f f:,(o ■ ■ ■ fm ■ fvM ■ ■ ■ f Pt (o ■ m) 

We conclude that in the leading order in f n , f ft , the calculation of the matrix elements 
of the matter operators is reduced to the substitution 

a n - e~ 2in tf n , n > 

a. n - e 2 **/* , n > 

with subsequent integration over £. This proves the relation between the matrix 
elements ( p3[) and their classical counterparts 



26 



Let us turn to the matrix elements ([35]). Since p(+) and afjj. commute with Om 
and with the DDF operators, the operator d + X^ reduces to 

d + x[ +) = ~7>« (75) 

when sandwitched as in eg . (|35"D . To find the matrix elements involving d + X <y L ' one 
notices that in the leadig order in f n , 

^ pi \L m O M \^p) = (76) 

because the commutator of L m and Om does not contain f n , (recall that \^p) 
and l^p') are physical states). Equation ( [75| ) and just established relation between 
the matrix elements (^) and classical correlators (|37|) immediately imply the desired 
relation between the matrix elements (|35D and their classical versions (|38"D . This 
relation can of course be obtained by an explicit calculation. 
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